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(No) dynamical constraints on the mass of the black hole in two ULXs 
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We present the preliminary results of two Gemini campaigns to constrain the mass of the black hole in an ultraluminous 
X-ray source (ULX) via optical spectroscopy. Pilot studies of the optical counterparts of a number of ULXs revealed two 
candidates for further detailed study, based on the presence of a broad He II 4686 A emission line. A sequence of 10 
long-slit spectra were obtained for each object, and the velocity shift of the ULX counterpart measured. Although radial 
velocity variations are observed, they are not sinusoidal, and no mass function is obtained. However, the broad He II line 
is highly variable on timescales shorter than a day. If associated with the reprocessing of X-rays in the accretion disc, its 
breadth implies that the disc must be close to face-on. 
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1 Introduction 

The key uncertainty responsible for driving the study of ul- 
traluminous X-ray sources (ULXs) over the past decade is 
the mass of the compact objects powering this extraordinary 
phenomenon. This unknown has been addressed by various 
ingenious methods, drawing evidence from across the elec- 
tromagnetic spectrum, many of which are discussed else- 
where in these proceedings. (A separate discussion of many 
of these methods, and their results, is presented in Zampieri 
& Roberts 2009). A consensus has emerged that ULXs are 
powered by accretion onto a black hole; but the question of 
the mass of these black holes remains unanswered. 

The reason that ULX masses remain mired in contro- 
versy is that none of the currently utilised methods pro- 
vides a direct, unambiguous measurement of the black hole 
mass. This is perhaps best exemplified by X-ray spectral 
analyses. A great deal of progress has been made in recent 
years on the basis of Chandra and XMM-Newton spectra, 
firstly identifying soft excesses consistent with the cool ac- 
cretion disc signature one would expect from a ~ 1000 A/© 
intermediate-mass black hole (e.g. Miller et al. 2003; Miller, 
Fabian & Miller 2004). Latterly this interpretation has been 
strongly challenged by the detection of a spectral break at 
energies of a few keV, identified in the best quality XMM- 
Newton data for a wide range of ULXs. This implies that 
ULXs are operating in an unfamiliar spectral state, most 
likely associated with super-Eddington processes (Stobbart, 
Roberts & Wilms 2006; Roberts 2007; Gladstone, Roberts 



& Done 2009; also Gladstone, these proceedings). This 'ul- 
traluminous state' appears to display the characteristic im- 
print of a strong outflowing wind, as predicted for super- 
Eddington emission (e.g. Begelman, King & Pringle 2006; 
Poutanen et al. 2007), and so implies that ULXs harbour 
small, < 100Af Q black holes. 

However, neither of the above examples provides a di- 
rect mass estimate and, worse still, for most ULX X-ray 
spectral data degeneracy is a problem as the quality is suf- 
ficiently poor that neither model can be rejectee^ Other in- 
direct methods suffer similarly - the QPOs detected in the 
power density spectra of NGC 5408 X- 1 can be used to in- 
fer the presence of an IMBH (e.g. Strohmayer & Mushotzky 
2009; also Strohmayer these proceedings), but this assumes 
both a specific type of QPOs and a sub-Eddington accre- 
tion state. Neither may be true for this source (Middleton et 
al. 2010; also these proceedings). Similarly, different mod- 
els for the optical colours and magnitudes of various ULX 
counterparts lead to a range of mass estimates (e.g. Cop- 
perwheat et al. 2006; Madhusudhan et al. 2009; Patruno & 
Zampieri 2010). It is clear therefore that we require a 'clean' 
test of ULX mass, untainted by model assumptions. 

The obvious way forward is to perform similar exper- 
iments to those that have a near four-decade heritage for 
Galactic systems: dynamical studies based on the co-orbital 
motion of the accreting compact object and its companion, 
donor star (see e.g. Charles & Coe 2006; Casares 2007). 
In such experiments one commonly obtains optical (and/or 
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1 Indeed, other spectral models may also be applied in some cases, 
for example slim disc spectra (e.g. Vierdayanti et al. 2006), or reflection- 
dominated spectra (Caballero-Garcia & Fabian 2010), to name but two. 
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UV/IR) spectra at a series of different epochs, and mea- 
sures the semi-velocity amplitude K and period P for the 
sinusoidal orbital motions, as traced out by shifts in the ob- 
served emission and/or absorption line wavelengths. Mea- 
surements are usually taken from periods when the donor 
star dominates the optical light, and used to infer a mass 
function f(M) that places a lower limit on the black hole 
mass, Mx- However, it is also possible to use emission fea- 
tures originating in the accretion disc to produce a mass 
function (e.g. Orosz et al. 1994; Soria et al. 1998) such that 



f{M) 



sin 3 i 



(M c + M x ) 2 
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2vrG : 



(1) 



where Mc is the mass of the companion star, Kx the 
semi-velocity amplitude of the black hole, and i the incli- 
nation of the orbital plane to our line-of-sight, thus placing 
limits on the mass of the black hole. 

We can take encouragement from the recent reports of 
mass functions for two extragalactic BHBs, M33 X-8 and 
IC 10 X-l (Orosz et al. 2007; Prestwich et al. 2007; Sil- 
verman & Filippenko 2008). However, most ULXs are at 
least three times more distant than these objects, with my ~ 
20.5 at best (Motch, these proceedings) and more typically 
my > 24 (Roberts, Levan & Goad 2008). Furthermore, 
many are located in complex fields, where their spectro- 
scopic signal could be confused with neighbouring nebu- 
losity and/or stars. Indeed, only a handful of ULXs might 
be accessible for mass function measurements with current 
facilities. Very few attempts have been made to date, and 
these have been unsuccessful. Kaaret & Corbel (2009) ob- 
tained 6 VLT/FORS observations of NGC 5408 X-l over 3 
days, but found no stellar absorption lines to base a mass 
function on. Pakull, Grise & Motch (2006) obtained multi- 
epoch data for NGC 1313 X-2 and did find an interesting 
velocity shift (Av ~ 380 km s _1 ) in the centroid of a broad 
He II line, but were unable to constrain a mass function from 
subsequent follow-up data (Grise et al. 2009). Hence the 
first mass function measurement for a ULX remains a tan- 
talising goal. Here, we detail the preliminary results of a 
new attempt to obtain the mass functions of two ULXs, Ho 
IX X-l and NGC 1313 X-2, using the Gemini observatory 
telescopes. 



2 Steps towards determining the dynamical 
mass of a ULX 

We have been developing a programme for the past few 
years, with the sole aim of obtaining a dynamical mass mea- 
surement for a ULX. The programme has three main steps: 
(i) identify the optical counterparts to ULXs on the basis of 
the most accurate available X-ray positions from Chandra, 
and HST imaging; (ii) obtain pilot optical spectroscopy to 
investigate the presence of useful spectral features; and (iii) 
undertake the radial velocity measurements campaign. 
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Fig. 1 Pilot optical long-slit spectra of the ULX coun- 
terparts. Data were taken by the GMOS instruments on 
Gemini-N (NGC 5204 X-l & Ho IX X-l) and Gemini-S 
(NGC 1 3 1 3 X-2), using the B600 gratings. The optical mag- 
nitudes of the counterparts (extinction-corrected Vegamags, 
in HST filters) and exposure times were: TO555 = 23.3, 3 hr 
(NGC 1313 X-2); m 606 = 22.3, 0.8 hr (NGC 5204 X-l); 
m 555 = 22.5, 1.5 hr (Ho IX X-l). 



2.1 Pilot spectroscopy results 

In the first step we surveyed nearby (d < 5 Mpc) ULXs 
with available Chandra and HST data, and selected rela- 
tively bright (my < 23.5) and isolated objects for further 
study. In Fig. [T] we show the pilot spectra, obtained using 
the GMOS instruments on the Gemini telescopes, for three 
X-ray luminous (Lx > 5 x 10 39 ergs -1 ) ULX counter- 
parts. Each has previously been identified in the literature, 
with positions shown by e.g. Liu, Bregman & Seitzer (2004, 
NGC 5204 X-l) and Ramsey et al. (2006, NGC 1313 X-2 
& Ho IX X-l). 

Fig. Q] shows the pilot spectra are dominated by a rel- 
atively featureless continuum, and emission lines from the 
bubble nebulae known to surround each of these three ob- 
jects (Pakull & Mirioni 2002). However, one interesting fea- 
ture is seen in two of the spectra: a broad He II 4686 A line 
is evident in the spectra of NGC 1313 X-2 and Ho IX X- 
1 . It is very plausible that this line could originate from the 
reprocessing of X-rays in the outer regions of the accretion 
disc, so could be used to trace radial velocity variations (as 
noted for NGC 1313 X-2 by Pakull et al. 2006) and there- 
fore place limits on the black hole mass using Eq.Q] 



2.2 New Gemini spectral monitoring campaigns 

Ten follow-up Gemini observations were obtained for each 
of NGC 1313 X-2 (2.5 hr per observation on Gemini-S) and 
Ho IX X- 1 ( 1 .5 hr per observation on Gemini-N) in semester 
2009B. Nine of the observations of NGC 1313 X-2 were 
performed over a 13-day period in December 2009, with a 
view to sampling over the known ~ 6 day photometric pe- 
riod of this ULX (Liu, Bregman & McClintock 2009). The 
observations of Ho IX X-l were split into two blocks of 
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Fig. 2 Observed radial velocity shifts of the ULX coun- 
terparts relative to the bubble nebulae, measured from the 
broad He II 4686 A line, over the observation campaign. 
Top panel: Ho IX X-l. Lower panel: NGC 1313 X-2. 



5 observations, in late December 2009 and February 2010. 
Details of the data analysis will be presented by Gladstone 
et al. (in prep.). In brief, spectra were extracted for both the 
counterpart, and for the surrounding nebula. The redshift of 
the region containing the ULX was constrained from the 
(off-ULX) nebular emission line spectrum for each obser- 
vation, and this was then used as a fiducial marker to search 
for relative changes in the He II 4686 A line centroid from 
the contemporaneous counterpart spectrum. 



3 Provisional results 

The results we present here are based on an initial analysis 
of the data, and focus primarily on NGC 1313 X-2. A fur- 
ther, more complete analysis is in preparation by Gladstone 
et al. and Roberts et al. 

The preliminary results show that the campaign was suc- 
cessful on one count: we detected the anticipated shifts in 
the He II 4686 A line with respect to its local vicinity (see 
Fig. O, measuring velocity shifts of ±100 km s _1 for Ho 
IX X-l, and up to - 200 km s" 1 for NGC 1313 X-2. How- 
ever, the data was not consistent with sinusoidal variations 
in either case. Simple sinusoid fitting to both datasets re- 
sulted in poor fits, with P = 1.73 days, K — 77 km s _1 
and xl ~ 5 for Ho IX X-l; and P = 3.01 days, K = 
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Fig. 3 He II 4686 A line profiles from the NGC 1 3 1 3 X-2 
data. The best fitting Gaussian profile is shown as a solid 
blue line, and the data residuals to this fit as a dashed blue 
line below the data. The vertical dashed and dotted lines 
show the rest wavelength and the line centroid, respectively. 
The relative flux is shown with the peak line flux normalised 
to unity for all good (> 9a) line detections. The last two 
detections have lower statistical significances, of 4 and 6cr 
respectively. 

61 km s- 1 and xl ~ 4 for NGC 1313 X-2. Given the lack 
of evidence for sinusoidal (i.e. orbital) variations in the data, 
we must therefore conclude that the data does not provide a 
new, solid constraint on the mass function for either object. 

Despite the lack of mass function measurements, the 
data does reveal interesting behaviour from the ULX coun- 
terpart. One notable feature is the highly variable nature 
of the broad He II 4686 A line. In Fig. [3] we show the 
profile of this line in the 10 separate observations of NGC 
1313 X-2. Over these observations, its FWHM varies from 
a peak of ~ 900 km s _1 down to the instrumental resolu- 
tion (< 290 km s _1 ), with variations of factor 2 seen in 24 
hours. This dramatic variability implies at the minimum that 
the broad He II emission must be originating within 24 light- 
hours of the ULX. 

Such variations might come from X-ray reprocessing in 
the outer regions of the accretion disc. If so, the broad lines 
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Fig. 4 Mass limits on the black hole in NGC 1313 X-2, assuming the He II 4686 A line originates in the accretion 
disc. We use the 2a upper limit of 182 km s _1 derived from the rms scatter of the velocity shift, companion stellar mass 
estimates based on Patruno & Zampieri (2010), and calculate lower limits on the mass (for a range of inclinations, as per 
the legend) based on Eq.Q] 



tell us that the disc cannot be perfectly face-on (or else we 
would see no line-of-sight velocity variations). This means 
that there should be some information on the black hole or- 
bit hidden within the radial velocity data: we therefore use 
the 2cr upper limit derived from the scatter in radial veloc- 
ity measurements for NGC 1313 X-2 to derive lower lim- 
its on its black hole mass for a range of orbital periods in 
Fig |U However, as Kaaret & Corbel pointed out for NGC 
5408 X-l, such small measured velocities are inconsistent 
with the expected velocities of material in the accretion disc 
C$> 1000 km s _1 ). A possible solution is that only a small 
component of the accretion disc rotational velocity is in the 
line-of-sight, i.e. the disc is very close to face-on (i ~ 0°). 
We see in Fig [4] that even at i — 30° there is little or no 
lower limit on the black hole mass, other than at very short 
period^]. Hence if the disc is close to face-on, little can be 
said about the black hole mass. 



4 Conclusion 

The initial analyses of our campaigns to determine mass 
functions for two ULXs betray no strong evidence for pe- 
riodic velocity variations from either object; and so no mass 
function is forthcoming. We do see some interesting phe- 
nomenology in the ULX counterparts, including strong vari- 
ability in the He II 4686 A line indicating it originates very 
close to the ULX. The FWHM of the broad line indicates 
that the ULX may be close to face-on; interestingly, both 
objects are X-ray luminous (Lx > 5 x 10 39 erg s _1 ), and 
models of super-Eddington discs predict that their observed 
flux will be highest close to this line-of-sight (Mineshige, 
these proceedings). This is perhaps more circumstantial ev- 
idence for models of ULXs as super-Eddington stellar-mass 
systems rather than IMBHs. In the meantime, we still await 
the first determination of the mass function for a ULX. 



2 Although removing the last two, low significance radial velocity 
points, significantly reduces the rms scatter to 44 km s _1 , and this can 
lead to some interesting constraints, e.g. for a 16 Mq companion and 30° 
inclination, Mbh > IOMq over all periods up to 10 days. 
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